Anomalous Quasiparticle Lifetime in Graphite: Band Structure Effects 
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We report ab initio calculation of quasiparticle lifetimes in 
graphite, as determined from the imaginary part of the self- 
energy operator within the GW aproximation. The inverse 
lifetime in the energy range from 0.5 to 3.5 eV above the 
Fermi level presents significant deviations from the quadratic 
behavior naively expected from Fermi liquid theory. The de- 
viations are explained in terms of the unique features of the 
band structure of this material. We also discuss the experi- 
mental results from different groups and make some predic- 
tions for future experiments. 
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The Fermi liquid theory ||l| (FLT) has been one of the 
most successful paradigms to describe the metallic state 
of condensed matter. According to it, the excitation 
spectrum in a metal can be understood in terms of the 
quasiparticle (QP) concept: Quasiparticles are in one-to- 
one correspondence with the one-electron excitations of 
a non-interacting Fermi system, but their properties are 
modified by the interactions. In a Fermi liquid, for ener- 
gies very near the Fermi level, quasiparticles have a finite 
lifetime (r) , which is inversely proportional to the square 
of the excitation energy ||l|,pi, i.e. t oc {E — Ep)~'^. It 
is tempting to claim a breakdown of the paradigm when 
deviation from this behavior is found for a metallic sys- 
tem. Recent experiments by Xu and coworkers B on 
graphite have reported an anomalous behavior of the QP 
lifetime for E — Ep in the range of 0.4 to 2 eV. These 
authors observed that t oc {E — Ep)~^. In this letter, 
we argue that, when the band structure of graphite is 
taken into account to compute the QP lifetime within 
FLT, important deviations from {E — Ep)~'^ behavior are 
found in the experimentally accessible energy range. We 
have found the QP lifetime to be strongly dependent on 
the wave vector direction and that, when averaged over 
the first BriouUin zone (BZ), the energy dependence of r 
cannot be fitted to a simple power law over the consid- 
ered energy range, from 0.5 to 2.5 eV. All these features 
find a simple explanation in terms of the dependence of 
the electron self-energy on the peculiarities of the band 



structure of graphite. 

In the experiments carried out by Xu et al. the sam- 
ple (in this case. Highly Oriented Pyrolytic Graphite, 
HOPG) is excited by a laser pulse and a second delayed 
pulse is used to probe the population of the excited states. 
Varying the delay time between the first (pump) and 
the second (probe) pulses provides a powerful, though 
still poorly understood, spectroscopic technique known 
as two-photon photoemission (2PPE). Subsequent mea- 
surements by Ertel and coworkers y, using very similar 
techniques, have partially confirmed these results. They 
found that r oc {E — Ep)~°' with a « 1.2 — 1.3 provided a 
better fit to their 2PPE data, for 0.2 eV < £■ - £^f < 1-6 
eV. It should also be pointed out that the absolute value 
of r seems to depend on the method chosen to fit the 
2PPE signal as well as the way in which the transport 
processes at the surface are accounted for Q. However, in 
both experiments the deviation from the {E — Ep)~'^ de- 
pendence (found in many simple and noble metals P,pl) 
seems clear. 

Graphite is a layered semimetal. The coupling between 
layers (i.e., along the c-dircction) is much weaker than 
within the layers (i.e. in the a6-plane), where strong 
bonds of covalent type keep the atoms arranged in a 
honeycomb lattice with two atoms per unit cell. This 
arrangement makes the valence and conduction 7r-like 
bands touch at the corners of the hexagonal first BZ. 
Interlayer coupling in an ABAB... stacking sequence, 
however, enhances this overlap and shifts the Fermi level, 
thus leading to the formation of electron and hole pock- 
ets responsible for the semimetallic character of graphite. 
Gonzalez et al. [[7| have put forward an explanation that 
relates the apparent linear dependence of the inverse QP 
lifetime to the structure of the degenerate 7r-bands near 
Ep and the weak coupling between layers. These two 
factors result in a small density of states at Ep. As 
consequence, the screening length in graphite is very 
large ^j^, a fact that strongly affects QP lifetimes. How- 
ever, the theory developed by Gonzalez and coworkers 
relies on a simplified model for the band structure of 
graphite, which assumes that the 7r-bands disperse lin- 



early with the wave vector parallel to the a6-plane (k||) 
up to an energy of w ±3 eV around Ep. Such a sim- 
plified band structure is obtained as a low energy limit 
of a tight-binding fit of the 7r-bands g| , and completely 
neglects interlayer hopping. Our and previous band- 
structure calculations, however, indicate that, along some 
k-directions, the dispersion is strongly non- linear already 
for E—Ep ~ 1-5 eV. Moreover, hopping between graphite 
layers, although smaller than intralayer hopping, is not 
negligible (w 0.25 eV) and leads to a splitting of the tt- 
bands at all energies. 

Motivated by these findings and the experiments men- 
tioned above, we have evaluated the QP lifetime account- 
ing for the band structure of graphite (lO| . The lifetime 
for a QP with wave vector k and band index n was ob- 
tained from the imaginary part of the electron self-energy 
nE) 1,0: 
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where we have neglected the renormalization of the QP 
wave function and energy by taking the expectation value 
over the corresponding LDA orbital |nk) and setting 
E — -Enk, the LDA eigenvalue. S(-E) was obtained within 
the GW approximation ||l2|,^. Since we are interested 
in a spectral property, namely the QP lifetime, it is suffi- 
cient to stop at this level [Q . Thus the screened electron- 
electron interaction can be calculated within the random 
phase approximation (RPA) [ [l5[ in terms of the dielec- 
tric matrix eGG'(q, w) [Hll|. Finally, since the 2PPE 
experiments can only probe the energy dependence in 
the relaxation time of the excitations generated by the 
pump pulse, we have averaged t~^ over the states with 
energy E = ii'nk to obtain an energy-dependent inverse 
lifetime, t-^{E) §. 
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FIG. 1. Graphite 7r-band structure near some relevant 
points M, L and K in the energy range considered in this 
work. 



The reason that the average QP lifetime cannot be 
fitted to a simple power law over the considered energy 
range is because of the smallness of r^^ in the neighbor- 
hood of the M and L points of the first BZ, where the 
TT-bands exhibit a saddle point (see Fig.l). This band 
topology forces a QP excited near these points to de- 
cay into a lower energy state by transferring a relatively 
large momentum (q||) parallel to the a6-plane. On the 
other hand, the large screening length referred to above 
favors small q|| electron-electron scattering and, at the 
same time, the band structure below 1.5 eV leads to a 
small density of low energy excitations available at large 
momenta. This point is illustrated in quantitative detail 
by Fig. 2, where it can be seen that the spectral weight 
of the energy loss function in the low frequency region 
rapidly decreases with increasing |q|||. This feature is 
also present in Imeoo(q, w) (not shown here for simplic- 
ity), which exhibits a similar structure, revealing that the 
low energy excitations are indeed electron-hole pairs and 
not plasmons. Thus the main decay channel for QP's is 
provided by electron-hole pair creation and not by plas- 
mon excitation as speculated in Ref. ||] (see, however, 
Refs. JzlJlSl). The coarseness of the grid used in our cal- 
culation precludes us from obtaining the lifetime below 
0.5 eV where we expect other decay channels such as 
phonons or low energy plasmons p^ to be very impor- 
tant as well. 
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FIG. 2. Energy loss function, Ime(,Q(q, 
for different values of q. 
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In Fig. 3 our results for r~J are plotted vs. the LDA 
energies, E^-k, for several k-directions: K ^ T^ and 
M ^ V. Notice that the value of the QP lifetime de- 
pends strongly on k. As discussed above, the unoccupied 
TT-bands possess a saddle point near M (cf. Fig.l), where 
the QP lifetime increases dramatically. This explains the 
rapid drop of r~J along M ^ F, for i? - Ep « 1.5 eV. 
Above this energy, the states around M become avail- 



able for a QP to decay and t~^ becomes appreciable. As 
E — Ep increases further, more bands come into play and 
thus the lifetime aquires two values near E — Ep w 2.3 
eV, which correspond to states with similar energies but 
different band index, n. The splitting of the lowest en- 
ergy bands along the K ^ T point also has the same 
effect. This splitting is due to the interlayer hopping, 
t± w 0.25 eV, which lifts the degeneracy of the 7r-bands 
near Ep. 
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FIG. 3. Inverse QP lifetime t'^ 
two representative directions K -^ 
(squares). 
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As to the K ^ r direction, one can regard the behavior 
of T~j. as controlled by two parameters: The excitation 
energy, e = E — Ep « vp |k|| |, where vp is the Fermi ve- 
locity, and the interlayer hopping t± . If we focus on the 
average behavior of the lifetime along this direction, given 
the small splitting of the bands, we can drop the band 
index n. Let T^^{e,t±) be the QP inverse lifetime. On 
dimensional grounds, T~^(e,tj_) = g^e f{t±/e, g), where 
g oc e'^/vp is the dimensionless constant "charge" that 
characterizes the strength of the electron-electron inter- 
action |^,g] and f{x,g) is a scaling function. If t_L = 
then the model of Ref. m would apply for the QP states 
along K ^ r. Therefore, r^^ ~ g^ e, which implies 
f{0,g) — const. On the other hand, for t± finite the 
QP lifetime is given by Fermi liquid theory [pi so that 
T~^ r-^ g^e^ as e — > 0, and therefore f{x ^ l,g) ^ 1/x. 
This means that we can expect a crossover from a regime 
where r~^ has an approximate linear behavior with e 
to the quadratic behavior expected from FLT. The lat- 
ter can be hardly seen in our results, which only extend 
down to e « 0.5 eV, but, since we rely on the QP picture 
of FLT for our calculations, it should be recovered suffi- 
ciently close to Ep. On the other hand, for e well above 
ij_, we have approximately r^^ oc e along iiT ^ P, as seen 
in Fig. 3. Thus our results for this direction agree qual- 



itatively with the linear dependence found by Gonzalez 
et al. y\ on the assumption that e « WF|k||| but using a 
different approach. 
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FIG. 4. Average inverse QP lifetime t~^{E) vs. QP energy 
E: present calculation (circles) and experimental results from 
Ref. [3] (squares) and Ref. [18] (triangles). 

Finally, in Fig. 4 we show the results for the inverse life- 
time of all the k-points considered and the resulting av- 
erage, which shows three different regimes: r^^ increases 
monotonically between 0.5 and 1.0 eV. From there and 
up to 1.5 eV, it remains approximately constant and even 
decreases due to the long QP lifetime found near M and 
L at E~Ep « 1.5 eV, which was explained above. Above 
1.5 eV r"^ increases rapidly and, as the energy increases, 
band structure effects should become less and less impor- 
tant. Also, for the sake of comparison we have plotted 
the experimental results reported in Ref. [p[, and new 
experimental results of Moos et al. [M . A word of cau- 
tion is in order as the precise interpretation of what is 
measured in these experiments in terms of QP lifetime 
is far from clear [H,H,n8|; results from different groups 
using very similar techniques and interpretations do not 
always coincide. However, it is worth commenting on 
the fact that the results by Xu and coworkers do not 
show any dip around 1.5 eV that is predicted by our 
results but instead increase monotonically in the whole 
energy range. At first, one could blame the approxima- 
tions made in the present calculation. Certainly, if the 
renormalization of the QP energy is not neglected it can 
be expected that the position of the saddle points near 
L and M may be slightly shifted upwards in energy and 
this would extend the range over which t^^{E) increases 
monotonically with E. Also, detailed photoemission ma- 
trix elements, not considered here, may affect the mag- 
nitudes of the decay rates. Other QP decay mechanisms 
such as electron-phonon coupling should be considered 
as well. However, the results of Moos et al. [nSl on pris- 



tine graphite surfaces show a dip in t~^ around 1.5 eV 
as predicted by the calculations reported here. Further- 
more, they also observe a strong dependence on the sam- 
ple preparation, which might explain discrepancies with 
the results from other groups. 

In conclusion, we have presented the results of a cal- 
culation of the lifetime, r, which takes into account the 
band structure of graphite. We find important deviations 
from the r oc (£' — Ep)^"^ naively expected from Fermi 
liquid theory. The deviations, however, find a natural 
explanation given the peculiarities of the band structure 
of graphite. The saddle point near the M point leads to 
a very large r (alternatively a very small t~^), whereas 
the linearly dispersing bands near along the iiT ^ F di- 
rection lead to r oc (£■ — Ep)~^ for energies well above 
the interlayer hopping w 0.25 eV, in agreement with pre- 
vious work . Finally, when averaged over the BriouUin 
zone, we predict that t^^{E) cannot be fitted by a simple 
power law as it exhibits a dip near 1.5 eV. This behavior 
has been observed in recent experiments |18| . 
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